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INTRODUCTION 


The  U.S.  Navy  has  initiated  a  significant  effort  to  develop  robust  modeling  tools  to  predict  structural 
response  of  high-strength,  high-toughness  ship  and  submarine  hull  steels  to  underwater  explosive  (UNDEX) 
loading.1  Naval  structures  subject  to  UNDEX  loading  experience  a  spectrum  of  “high-rate”  loading  conditions, 
encompassing  dynamic  buckling  at  small  strains  to  large  strain  deformation  prior  to  hull  rupture.  Both  loading 
conditions  realize  strain  rates  prior  to  “localization”  and  failure,  which  are  on  the  order  of  100  s”1.  The  fidelity  of 
structural  response  prediction  relies  heavily  on  the  accuracy  of  the  constitutive  model  and  supporting  calibration 
data  sets.  In  particular,  constitutive  strength  models  must  accurately  describe  yielding  and  small  strain  behavior, 
in  which  a  rapidly  evolving  tangent  modulus  highly  influences  the  onset  of  dynamic  buckling,  to  large  strain 
bulging  and  plastic  deformation  where  stress  states  exist  between  uniaxial  and  moderate  stress  triaxialities.  In  the 
latter  case,  failure  may  initiate  at  large  plastic  strain  levels,  such  as  those  characterized  by  zpi— >  1.0,  at  which  the 
strain-hardening  behavior  must  be  accurately  represented.  As  such,  the  U.S.  Navy  is  currently  interested  in 
investigating  and  quantifying  the  constitutive  stress-strain  response  of  ship  and  submarine  hull  steels  (HY-80, 
HY-100  and  HSLA-100)  and  their  associated  weld  systems. 

The  constitutive  stress-strain  response  of  a  U.S.  Navy  HY-100  steel  is  characterized  employing  a  three-term 
mechanical  threshold  stress  (MTS)  strength  model  approach.  MTS  is  a  semi-empirical  description  of  the  plastic 
deformation  (stress-strain)  behavior  of  (primarily)  metals  and  metallic  alloys,  and  is  developed  on  the  basis  of  two 
principal  concepts.  Flow  stress,  o,  is  a  metric  of  the  “microstructural”  state  of  a  material,  i.e.,  stress  is  a  state 
variable.  A  mechanical  threshold  stress,  a,  is  also  defined  as  the  flow  stress  at  a  temperature  of  absolute  zero 
(i.e.,  0  K)  such  that  deformation  can  be  quantified  in  the  absence  of  thermally  activated  processes.  To 
accommodate  the  effects  of  thermal  work,  the  mechanical  threshold  stress  is  scaled  based  on  an  Arrhenius 
description  of  thermally  activated  deformation.  For  brevity,  only  the  salient  features  of  the  MTS  approach  will  be 
discussed.  The  interested  reader,  however,  is  directed  to  the  background  and  discussions  of  the  approach  by 
Kocks  et  al.2'5 

The  mechanical  threshold  stress,  6,  is  commonly  expressed  as  a  superposition  of  flow  stress  contributions, 
each  component  of  which  characterizes  a  relevant  deformation  mechanism: 

-  =  ^-+S^-  (i) 

^  JU  \i 

where  p  is  the  temperature-dependent  shear  modulus  and  a  represents  the  mechanical  threshold  stress  for 
athermal,  or  temperature-independent,  processes  (subscript  a)  and  for  the  jth  individual  thermally  activated 
deformation  process  (subscript  j).  Phenomenologically,  athermal  deformation  processes  primarily  consider  the 
“Hall-Petch”-like  grain  boundary  strengthening  contribution,  whereas  the  individual  thermally  activated  processes 
encompass,  for  example,  dislocation  interactions  with  solute  (interstitial  and  solid  solution)  atoms  and  second- 
phase  precipitates  and  dislocation  interactions  with  other  dislocations.  Recall  the  mechanical  threshold  stresses, 
6a  and  &  j ,  are  the  flow  stress  components  from  athermal  and  thermally  activated  processes  in  the  absence  of 

thermal  work.  The  strength  model  incorporating  the  MTS  concept  and  rate/thermal  effects  is  outlined  below. 
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Temperature  (thermal)  effects  are  incorporated  into  the  MTS  strength  model  framework  by  a  temperature 
(and  strain-rate)  dependent  scaling  factor,  S,  where: 


(2) 


where  p  and  go  are  temperature-dependent  and  0-K  shear  moduli,  respectively,  and  S  is  the  ratio  of  temperature- 
dependent  “flow”  stress  to  the  mechanical  threshold  stress  \pj/Oj}-  The  scaling  factor,  5,  is  derived  from  an 

Arrhenius  expression  relating  strain  rate  to  activation  energy  and  temperature: 


e  =  80Jexp  J/kT 

V  J 


assuming  a  stress-dependent  activation  (free)  energy,  AG: 


AG  j  =  AG| 


=  80,j^b: 


1- 


(  n  : 


\PJ 


'a 


V 


(3) 


(4) 


where  g0  is  a  normalized  activation  energy,  b  is  the  burgers  vector,  p  and  q  are  empirical  constants  related  to  an 
“obstacle  profile,”  eis  strain-rate,  e0  is  a  reference  strain-rate,  k  is  the  Boltzmann  constant,  and  T  is  absolute 
temperature.  Equations  (3)  and  (4)  can  be  manipulated  and  re-expressed  such  that: 


Sj  = 


i- 


kT 


SojMb 


-In 


L  V*. 


1  /pj 


'0,7 


(5) 


where  S  is  bounded  and  obeys  0  <  S  <  1 . 

The  role  of  temperature  (thermal  energy)  on  flow  stress  is  effectively  illustrated  by  equations  (2)  and  (5). 
Increasing  temperature  reduces  the  necessary  mechanical  energy  (stress)  for  plastic  deformation.  In  addition,  the 
equations  comprising  (5)  illustrate  that  the  temperature  dependence  of  flow  stress  also  gives  rise  to  explicit 
(strain)  rate  sensitivity,  viz.,  equation  (3).  In  other  words,  flow  stress,  o,  is  not  a  unique  quantity  with  respect  to 
either  strain-rate,  e ,  or  temperature,  T.  Therefore,  numerous  combinations  of  strain  rate  and  temperature,  of 
which  many  are  experimentally  tractable,  result  in  an  identical  value  of  flow  stress.  A  practical  application  of  this 
observation  is  the  potential  to  simulate  high  strain  rates/ambient  temperature  tests  via  low  strain  rates/low 
temperature  tests. 

Plastic  strain  is  not  explicitly  defined  within  the  MTS  strength  model  framework,  as  per  equation  (1). 
Plastic  strain  is  characterized,  implicitly,  in  terms  of  thermally  activated  dislocation-dislocation  interactions  and, 
specifically,  through  the  mechanical  threshold  stress,  a8  .  This  approach  to  interpreting  plastic  strain  arises  from  a 
viewpoint  that  the  evolution  of  the  dislocation  density  between,  for  example,  two,  constant  structure  regimes — 
characterized  by  ae  \  and  ae  2 — is  described  by  plastic  strain.  This  description  allows  changes  in  loading  path  to 
be  accurately  represented,  e.g.,  strain  rate  jumps  in  which  the  flow  stress  changes  from  )  to  a(e2)  where 
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e,  *e2 .  A  phenomenological  background  and  method  to  determine  plastic  strain  over  an  increment  of  6e  is 
described  below. 

The  relationship  between  plastic  strain  and  structure  evolution  (dislocation  density)  is  probably  best 
represented,  in  general,  by  considering  dislocations-dislocation  interactions  involving  the  entanglement  of 
dislocations  (storage/accumulation)  and  their  mutual  annihilation  (cross-slip/dynamic  recovery).  Microscopically, 
these  interactions  are  manifested  by  the  macroscopic  strain-hardening  rate,  0. 

6  =  dd/e  =  dh~dr  (6) 

where  Qh  and  0r  are  correlated  to  contributions  to  strain-hardening  from  dislocation  accumulation  and  dynamic 
recovery  respectively.  Note  that  dynamic  recovery  and  dislocation  accumulation  are  competitive  processes 
leading  to  “strain-softening”  and  “strain-hardening”  respectively. 

Physically  based  descriptions  for  0/,  and  0r  exist;  however,  the  strain-hardening  rate  is  frequently  described 
empirically  by: 


0  =  dd%i8  =  eO  IWM 


(7) 


where  0O  refers  to  the  stage  II  strain-hardening  rate  and  f(de )  characterizes  an  empirically  derived  dynamic 
recovery  rate.  Stage  II  hardening  is  typically  presumed  as  fairly  insensitive  to  temperature  and  rate.  At  small 
deformation  levels,  the  strain-hardening  rate  is  approximately  the  stage  II  hardening  rate,  i.e.,  0  ~  0O,  and  the 
dislocation  structure  is  characterized  by  a  small  value  of  de.  With  increasing  deformation  the  contribution  of 
may  be  substantial  such  that  f(ae)  1.  In  the  limit  of  f(6e)  =  1,  the  strain-hardening  rate  is  zero  and  a 
condition  of  constant  flow  stress  is  achieved,  i.e.,  flow  stress  saturation.  In  short,  the  saturation  stress  response  can 
be  interpreted  as  a  balance  being  achieved  between  the  rates  of  accumulation  and  recovery;  the  net  dislocation 
population  (dislocation  density)  neither  increases  nor  decreases  with  additional  deformation. 

Equation  (7)  can  be  rewritten  and  integrated  over  an  appropriate  ae  interval  to  calculate  plastic  strain.  The 
specific  form  of  s  is  dependent  on  the  expression  describing  the  strain-hardening  rate.  The  strain-hardening 
response  of  HY-100  (and  other  alloys)  is  given  by  a  modified  Voce  empirical  hardening  rule— the  tank  rule.  The 
tank  rule  captures  the  flow  stress  saturation  behavior  observed  in  OFE  Cu,  and  has  been  similarly  employed  to 
represent  the  strain-hardening  behavior  of  HY-100  (HY-80)  steel.  The  tank  rule  is  given  below: 


e=— -6ol 

dt 


tanh 


1-- 


a- 


Oe 


des(e,T)J 


tanh(a) 


(8) 


where  o  is  an  empirical  “best-fit”  constant  which  dictates  the  rate  at  which  saturation  is  achieved,  acs  is  a 
temperature-  and  rate-sensitive  saturation  stress,  represents  a  flow  stress  contribution  from  dislocation 
accumulation/annihilation,  and  0O  is  an  experimentally  determined  stage  II  strain-hardening  rate.  The  saturation 
stress,  Oc$ ,  is  derived  from  the  saturation  threshold  stress,  (TcSo ,  by  : 
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where  the  athermal  component  oa=aa,  and  subscripts  i  and  e  refer  to  an  “intrinsic”  barrier  to  thermally 
activated  dislocation  motion  and  dislocation-dislocation  interactions  and  microstructural  evolution  with  increasing 
deformation,  respectively.  Athermal  contributions  to  flow  stress  can  be  associated  with  dislocation  interactions 
with  martensite/bainite  lath/plate  boundaries,  i.e.,  Hall-Petch-type  behavior.  The  intrinsic  barrier  description 
encompasses  numerous  thermally  activated  deformation  processes.  Such  processes  are  of  principal  interest  in 
steels  and  may  potentially  include  such  mechanisms  as  overcoming  the  Peierls  barrier,  dislocation  interactions 
with  vacancies,  dislocation  interactions  with  solute  (interstitial  and  solid  solution)  atoms,  viz.,  carbon,  oxygen, 
and  nitrogen,  and  potentially  dislocation-carbide  (dislocation-particle  interactions).  Dislocation  interactions  with 
other  dislocations  are  also  thermally  activated.  However,  unlike  the  static  (constant  structure)  description  of 
microstructure  given  by  the  intrinsic  barrier,  e.g.,  “yield”  behavior,  the  Oe  term  describes  the  dynamic  evolution 
of  the  microstructure  with  increasing  deformation,  i.e.,  strain-hardening  behavior. 
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EXPERIMENTAL 


The  HY-100  steel  investigated  was  received  in  the  form  of  25-mm-thick  plate.  HY-100  is  a  quenched  and 
tempered  NiCrMo  steel.  The  nominal  plate  composition  is  0.16C-2.62Ni-1.32Cr-0.25Mo-0.26Mn-0.14Cu-0.22Si- 
0.008P-0.009S.  The  plate  was  austenitized  at  1,195  K  for  60  minutes  followed  by  a  water  quench.  Subsequently 
the  plate  was  tempered  at  910  K  for  90  minutes  followed  by  water  quenching.  The  HY-100  steel  plate  was 
received  and  tested  in  this  heat-treated  condition. 

The  MTS  strength  model  parameters  were  determined  based  on  uniaxial  compression  test  results  from 
HY-100  steel.  Solid  cylindrical  test  specimens  were  machined  from  as-received  plate  in  the  short-transverse 
orientation  (Figure  1).  Samples  were  6.35  mm  in  diameter  and  had  a  1:1  height-to-diameter  aspect  ratio.  Prior  to 
testing,  specimen  loading  faces  were  lubricated  with  molybdenum  disulfide.  Testing  was  performed  over  test 
conditions  encompassing  temperatures  and  strain  rates  between  77  and  298  K  and  0.001  and  ~2000  s  1 
respectively.  Multiple  reloads  were  performed  under  room-temperature,  quasi-static  loading  conditions  to  achieve 
plastic  strains  in  excess  of  1.0.  Upon  unloading,  typically  in  Ae  =  0.25  increments,  the  sample  was  relubricated  to 
minimize  frictional  effects;  a  minimal  degree  of  barreling  was  observed.  At  c  =  0.50  and  1.0  the  deformed 
specimen  was  remachined  to  maintain  a  1:1  aspect  ratio. 


a.  Uniaxial  compression  b.  Thin-wall  torsion 

Figure  1.  Uniaxial  Compression  and  Thin-Walled  Torsion  Specimens 


“Elevated  temperature”  testing  was  also  performed  at  dynamic  strain-rates  at  temperatures  of  473,  673,  and 
873  K.  Compression  specimens  were  heated  in  approximately  5  to  10  minutes  under  a  nominal  vacuum  of  10  3 
torr  to  the  test  temperature  and  thereafter  “soaked”  at  the  test  temperature  for  an  additional  5  minutes  prior  to 
testing.  However,  it  should  be  noted  that  o-e  data  acquired  at  these  elevated  temperatures  should  be  considered 
with  caution  as  the  additional  tempering  during  heating  and  soaking  may  result  in  a  minor  degree  of  softening  of 
the  original  HY-100  steel  martensitic/bainitic  microstructure. 
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Torsion  tests  of  thin-walled  tubes  at  298  K  and  1CF2  s_1  probe  the  large-strain,  shear-stress,  shear-strain 
response.  Specimens  were  machined  from  as-received  plate  such  that  the  torsion  axis  coincided  with  the 
transverse  plate  orientation  (Figure  1).  Shear  stress,  x,  and  shear  strain,  y,  for  thin-walled  tube  are  given  by  x  = 
T/2na2t  and  y  =  aQ/L,  respectively,  where  T  is  torque,  a  is  the  mean  radius  of  the  tube,  t  is  the  wall  thickness,  0  is 
the  angle  of  twist,  and  L  is  the  gage  length.  A  principal  motivation  for  conducting  torsion  tests  is  to  avoid  plastic 
instability,  as  occurring  in  uniaxial  tension  tests,  and  frictional  effects,  as  possible  in  uniaxial  compression  testing, 
prior  to  attaining  large  plastic  strains  (e  >  1 .0). 

In  addition,  texture  of  the  HY-100  steel  plate  was  probed  to  infer  the  extent  of  anisotropy  of  the  yield 
surface.  Specimens  were  sectioned  parallel  to  the  plate  surface  (L-LT)  in  the  vicinity  of  the  plate  surface  and 
center,  and  texture  was  determined  in  the  rolling  plane,  i.e.,  perpendicular  to  the  short-transverse  orientation 
(Figure  2).  The  texture  data  are  discretized  by  calculating  orientation  distribution  functions  (ODF),  using  the 
algorithm  of  Williams,  Imhof,  Matthies,  and  Vinel  (WIMV  analysis).6  This  ODF  information  subsequently 
motivates  a  yield  surface  determination  via  a  Taylor  and  Bishop-Hill  analysis.7  In  short,  this  procedure  calculates 
a  discrete  set  of  stress  values  comprising  the  yield  surface.  These  stresses  may  subsequently  be  fitted  to  an 
existing  yield  surface  formulation.  The  current  study  fits  the  discrete  stress  values  to  a  quadratic  “Hill-48”  yield 
surface  formulation,8 


2 


-ct33)2  +c(a33  _CTll)2  +  #(CT11  -  CT22  )2  +2Lct?3  +  2Ma?3  +2 Nof, 


■  a2  =  0 


(12) 


using  the  method  outlined  by  Maudlin  et  al.7  Expressed  in  terms  of  deviatoric  stress,  equation  (12)  can  be  re¬ 
expressed  as: 


-[(G  +  H)Sn  +  (F  +  h)S2 22  +{F  +  G)S332  -  2HSuS22  - 2 GSnSn  (13) 

2 

-  2FS22Si3  +  2LS23  +  2 MS2  +  2 AS2  ]  -  a2  =  0 


and  the  Hill  coefficients  can  be  expressed  in  terms  of  the  deviatoric  stresses  as: 


1 

S\2  1 

G  +  H 

1  1 

Si3  ‘ 

H  +  F 

S33  O' 

1 

S22 

F  +  G 

£ 
1  2  L 


(14) 
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RESULTS 


Yield  Surface  Determination 

The  yield  behavior  of  the  HY-100  steel  plate  is  inferred  as  isotropic.  The  coefficients  for  the  quadratic 
“Hill-48”  yield  function,  deduced  from  x-ray  crystallography  measurements  of  the  plate  material  and  the  Taylor- 
Bishop-Hill  analysis,  are  in  good  agreement  with  those  for  a  randomly  textured,  isotropic  material  deforming  by 
pencil  glide  (Table  I).  Pole  figures  and  an  inverse  pole  figure  measured  at  the  plate  center,  shown  in  Figure  3, 
indicate  m.r.d.  (multiple  of  random  distribution)  values  near  1.0,  i.e.,  a  near-random  texture.  Taylor-Bishop-Hill 
analysis  results  are  equivalently  represented  through  a  Tt-plane  representation,  shown  in  Figure  4.  The  yield 

J3 

surface  represented  in  Figure  4  is  symmetric  about  — (S 22  —5 n)  =  0 .  The  shear  sub-space  corresponding  to  the 

2 

HY-100  plate  material  is  spherical,  with  V3S23  =  V3.S,i2  =  V35^3 ,  because  L,  M,  and  N  are  essentially  identical. 


Table  I.  Quadratic  Yield  Surface  Coefficients 
for  a  Randomly  Textured  Isotropic  Material  and  the  Current  HY-100  Plate 


Texture 

F 

G 

H 

L 

M 

N 

Taylor 

factor 

Random  (pencil  glide) 

1.00 

1.00 

1.10 

3.24 

3.20 

3.21 

2.73 

Plate  (measured) 

1.00 

0.998 

1.13 

3.30 

3.25 

3.31 

2.73 

a.  Pole  Figure 

Figure  3.  Pole  Figures 


b.  Inverse  Pole  Figure 

of  HY-100  Steel  Plate  Measured  by  X-Ray 
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^3/ 2  (S22  —  S„) 


Figure  4.  Pi-Plane  Representation  of  the  Yield  Surface  for  As-Received 

HY-100  Steel  Plate  Material 

[Calculated  plate  material  surface  is  based  on  x-ray  texture  measurements. 
Randomly  textured  HY-100  plate  is  simulated  using  the  Taylor-Bishop-Hill  analysis  approach.] 


Compression  Test  Results 

The  MTS  strength  model,  equation  (11),  is  applied  to  quantify  the  experimentally  determine  mechanical 
response  of  HY-100  steel  shown  in  Figures  5  through  7.  Quasi-static  (10'3  and  10'1  s"1)  compressive  stress-strain 
responses  at  77,  188,  and  298  K  are  shown  in  Figure  5.  A  room-temperature,  dynamic  strain-rate  (e  >  103  s'1)  test 
result  is  given  in  Figure  6.  Strength  model  parameters  are  determined  from  the  data  based  on  a  procedure 
described  by  Chen  and  Gray.9  Albeit  the  MTS  strength  model  is  a  physically  based  constitutive  model  describing 
deformation  via  the  kinetics  of  dislocation  motion,  the  necessary  “material”  parameters  are  often  empirically 
derived. 

The  mechanical  data  presented  in  Figure  5  are  characteristic  of  the  deformation  behavior  of  body  centered 
cubic  (bcc)  metals  and  alloys.  An  appreciable  increase  in  the  0.002  offset  yield  strength  is  typically  observed  in 
bcc  metals  and  alloys  with  decreasing  temperature.  Specifically,  a  rapid  increase  in  the  quasi-static  yield  strength 
of  HY-100  steel  from  ~  800  to  -  1,300  MPa,  roughly  a  40%  increase,  is  recorded  between  298  and  77  K 
respectively.  This  behavior  is  commonly  attributed  to  the  resistance  to  dislocation  motion  afforded  by  the  Peierls 
barrier  and  interstitial  solute  atoms.  The  effectiveness  of  these  obstacles  is  very  temperature  (thermal  energy) 
sensitive.  The  marked  drop  in  yield  strength  between  77  and  188  K  illustrates  that  obstacles  are  increasingly 
ineffectual  in  impeding  dislocation  motion  with  increasing  temperature. 
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Figure  5.  Quasi-Static  Strain  Rate  Stress-Strain  Responses 
of  HY-100  Steel  between  77  and  298  K 


Figure  6.  Dynamic  (high)  Strain  Rate  Stress-Strain  Response 
of  HY-100  Steel  at  298  K 


IHTR2168 


Figure  7.  Dynamic  (high)  Strain  Rate  Stress-Strain  Response 
of  HY-100  Steel  between  473  and  873  K 


The  quasi-static  strain-hardening  rate,  also  manifested  in  Figure  5,  increases  by  a  factor  of  three  between 
298  and  77  K.  Specifically,  the  hardening  rate,  (0  =  do/ds),  increases  from  0298k  ~  600  MPa/ strain  to  077K  ~  2,200 
MPa/strain,  at  a  plastic  strain  of  0.10.  These  differences  in  strain-hardening  response,  although  not  completely 
understood,  are  typically  rationalized  on  the  basis  of  a  competition  between  dislocation  accumulation  (hardening) 
and  dislocation  annihilation  (softening).  Dislocation  accumulation  is  principally  motivated  by  strain  (or 
deformation)  considerations.  The  latter  (softening)  process,  being  sensitive  to  temperature,  for  example,  through 
dynamic  recovery  or  cross  slip,  is  less  effective  at  lower  temperatures  at  which  thermal  energy  and  the  driving 
force  for  diffusion  are  reduced.  Although  this  argument  provides  rationale  for  larger  strain-hardening  rates  at 
lower  temperatures,  i.e.,  77  versus  298  K,  the  lack  of  an  adequate  metallurgical-motivated  mathematical 
description  of  strain  hardening  precludes  an  accurate  quantification  of  this  process.  As  such,  the  MTS  strength 
model  (as  well  as  many  other  models)  relies  on  a  phenomenological  description,  viz.,  a  Voce-like  formulation,10 
of  strain  hardening. 

Dynamic  material  response,  as  illustrated  in  Figure  6,  is  evaluated  similarly  to  that  of  the  quasi-static 
behavior.  Due  to  the  lack  of  stress  equilibrium  at  small  strains  during  a  split  Hopkinson  pressure  bar  test,  flow 
stresses  are  typically  evaluated  at  plastic  strains  in  excess  of  0.01.  The  room  temperature  flow  stress  at  epl  ~  0.04, 
corresponding  to  -2,500  s"1  strain  rate,  is  approximately  950  MPa.  Evaluating  this  flow  stress  with  respect  to  the 
corresponding  quasi-static,  room  temperature  flow  stress  at  ep;  ~  0.04  (i.e.,  o0.04  ~  860  MPa)  results  in  an  estimated 
strain  rate  sensitivity,  m  =  A(ln  a)/ A(ln  s )  -  0.007 ,  a  typical  “m-value”  of  high  strength  steels.  At  larger  strain 
levels  of  £pi  -  0.10  and  0.15  the  m-values  are  essentially  identical,  viz.,  m  -  0.007,  to  that  at  epi  -  0.04.  A  strain¬ 
hardening  rate  is  roughly  estimated  as  0  ~  840  MPa/strain  at  a  plastic  strain  of  0. 10. 
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The  elevated  temperature  o-s  response  of  HY-100  steel  was  also  evaluated  under  dynamic  strain  rate 
conditions,  as  illustrated  in  Figure  7.  Reiterating,  the  mechanical  response  at  these  elevated  temperatures  is 
probably  not  strictly  that  of  HY-100  steel  as  a  result  of  the  limited  tempering  and  microstructural  (phase) 
evolution  sustained  during  heat-up  and  soaking.  Flow  stress  at  a  ep,  -  0.04  decreases  with  increasing  temperature 
from  -880  MPa  at  473  K  and  2,500  s'1  to  -580  MPa  at  873  K  and  3100  s"1.  The  hardening  rate  at  ep,  -  0.10 
exhibits  an  even  more  precipitous  drop  between  473  and  873  K  from  0  ~  788  MPa/strain  at  the  former 
temperature,  to  0  -  225  MPa/strain  in  the  latter. 


Torsion  Test  Results 


Room-temperature,  quasi-static  (y  -  10‘2  s'1)  torsional  stress-strain  response  of  HY-100  steel,  illustrated  in 
Figure  8,  exhibits  good  reproducibility.  Torsion  test  data  are  related  to  uniaxial  compression  data  by  representing 
both  data  sets  in  terms  of  effective  stresses  and  strains,  employing  the  second  invariant  of  the  stress  and  strain 
tensor.  These  second  invariants  can  be  expressed  as: 


and 


[(°1  -  °2  )2  +  (°2  -  )2  +  (°3  -  °1  )2  Y2 

£  =  -yr[(£l  -e2)2  +  (e2  ~e3)2  +(e3  -El)2. 


_  1 

CT=v^ 


(15) 


Figure  8.  Large  Strain  Stress-Strain  Response 
of  HY-100  Steel  Tested  in  Torsion  at  298  K  and  y  *  0.01  s'1 
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The  relationships  between  effective  stress  and  strain  and  uniaxial  compression  (or  tension)  and  torsion  are 
listed  below: 


Effective 

Uniaxial 

Torsion 

Quantities 

a 

o 

V3x 

e 

E 

y/v/3 

Effective  stress  and  strain  values  are  subsequently  “adjusted”  to  account  for  texture  by  comparing  yield 
strength  ratios  in  compression  and  shear,  033:0,3,  between  an  isotropic  von  Mises  material  and  the  current  HY-100 
plate.  Yield  strength  ratios  are  calculated  using  Hill  coefficients,  cf.,  equation  (14).  A  von  Mises  material  is 
characterized  by  F  =  G  =  H  =  \  and  L  =  M  =  N  =  3.  For  example,  the  yield  strength  ratio  between  compression 
and  shear,  033:013,  can  be  written  as: 


<733  _  2M 
O13  V  P  +  G 


(16) 


For  a  von  Mises  material,  this  ratio  is  V3  =  1.732.  The  ratio  exhibited  by  the  current  plate  material  is  1.804.  On 
this  basis  the  plate  material  is  expected  to  yield  at  a  stress  level  approximately  4%  higher  than  that  presumed 
assuming  a  von  Mises  yield  criterion.  Therefore,  the  0.002  offset  yield  strength  is  predicted  as: 

ct  =  1.04x-\/3xt  (17) 


For  the  current  plate  material  this  value  is  a  «  800  MPa. 


MTS  Parameters 

The  MTS  strength  model  parameters  are  determined,  based  on  a  procedure  described  by  Chen  and  Gray9, 
using  mechanical  data  from  Figures  5  through  7.  Evaluation  of  the  dynamic  material  response,  i.e.,  Figures  6 
and  7,  also  involved  considering  adiabatic  heating  effects.  Although  the  MTS  strength  model  is  a  physically  based 
constitutive  model  describing  deformation  via  the  kinetics  of  dislocation  motion,  the  necessary  “material” 
parameters  are  empirically  derived.  An  approach  to  determine  MTS  strength  model  parameters  is  given  below, 
with  the  corresponding  equations  and  parameters  for  HY-100  steel  summarized  in  Appendix  A. 

The  temperature-dependent  shear  modulus  for  HY-100  steel  is  given  by  Chen  and  Gray11  as  an  empirical 
relation  in  the  form: 


(18) 
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where  po,  D,  and  T0  are  in  essence  empirical  parameters.  The  values  of  p<>,  D,  and  T0  are  71,460  MPa,  2,910  MPa, 
and  204  K  respectively. 

The  athermal  stress,  o„,  contribution  in  equation  (11)  is  estimated  as  oa  -  40  MPa;  no  rigorous  analysis  was 
performed  to  arrive  at  this  value  of  oa.  However,  assuming  Hall-Petch-type  boundary  strengthening  as  the 
principal  active  athermal  strengthening  mechanism,  an  athermal  stress  on  the  order  of  40  MPa  is  not 
unreasonable.  For  example,  Tomita12  determined  for  a  0.2C-NiCrMo  steel  a  Hall-Petch  constant  characterizing 
the  martensite  packet  size  of  ky  =  0.578  MN  m"5'2.  Assuming  packet  sizes  on  the  order  of  100  pm,  the  Hall-Petch 
effect  is  k/T 1/2  ~  50  MPa. 


Thermally  activated  deformation  of  HY-100  steel  is  generically  characterized  within  the  context  of  the 
“intrinsic”  term,  <r(,  and  is  assumed  to  be  primarily  the  Peierls  barrier,  carbon  interstitial  solute  atom 
strengthening,  and  dislocation  strengthening  in  the  quenched  and  tempered  martensite.  A  normalized  activation 
energy  of  g0i  ~  1.161  and  a  mechanical  threshold  stress  of  a,  ~  1341  MPa  are  determinable  based  on  a  “Fisher” 
diagram13  (Figure  9  and  Table  II).  Specifically,  g0i  and  &t,  are  determined  when  mechanical  data  are  presented  in 
a  Fisher  plot  format: 


Pi 

f  ~  \ 

Pi 

Pi 

r  i  ] 

A 

\kT  In 

Ml 

J 

W 

W 

Joi  J 

| lb3 

l  *  II 

A 


(19) 


This  procedure  does,  however,  involve  making  an  a  priori  assumption  as  to  values  of  the  glide  obstacle  profile 
parameters,  p,  and  qh  and  the  reference  strain-rate,  eQi .  Empirical  observations  and  guidance  given  by  reference 
(2)  suggest  pi  and  q,  be  equated  to  1/2  and  3/2,  respectively,  whereas  eQi  is  chosen  as  1013  s"1.  A  single  set  of 
constant  values  of  p„  q„  and  e0j  for  HY-100  presumes  all  active  deformation  mechanisms  can  be  effectively 
represented  by  a  single  thermally  activated  deformation  mechanism  over  a  range  of  strain  rates  and  temperatures. 


Figure  9.  Modified  Arrhenius  (Fisher)  Plot  for  “Intrinsic”  Strengthening 
Component  of  the  MTS  Strength  Model  for  HY-100  Steel 

[Mechanical  threshold  stress,  di ,  and  normalized  activation  energy  value, ,  g0i,  are  determinable  from  this  plot.] 
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Table  II.  Parameters  Used  to  Calculate  gol  and  <r, . 


kl/nb3  ln(e 0/  / e) 

[(ay  -aa  )/n]P' 

T(K) 

e  (s’1) 

oy  (MPa) 

M  (MPa) 

0.10907963 

0.132357801 

77 

0.001 

1,288 

71,239 

0.200156472 

0.109237414 

188 

0.001 

875 

69,975 

0.252494489 

0.104077917 

298 

0.1 

782 

68,499 

0.276002618 

0.102736741 

298 

0.001 

763 

68,499 

0.196370459 

0.109411649 

298 

2500 

860 

68,499 

0.272794009 

0.099951907 

473 

2800 

700 

66,063 

0.3545421 

0.0924051 

673 

3000 

580 

63,240 

0.434362062 

0.089142151 

873 

3100 

520 

60,404 

The  structure  evolution  parameters,  comprising  equations  (2),  (5),  (8)  and  (9),  also  describe  thermally 
activated  deformation,  and  as  such  are  similarly  determined  by  analyzing  the  so-called  “Fisher”-plot  behavior. 
Inherent  complexity  and  a  lack  of  a  specific  understanding  in  the  evolution  of  the  dislocation  population  in  highly 
dislocated  martensitic/bainitic  steels  also  necessitates  relying  on  previous  experience  in  “simpler”  material 
systems.2,14  Therefore,  the  structure  evolution  glide  obstacle  profile  parameters,  pe  and  qE,  and  the  reference  strain 
rates,  e0c  and  e0eS ,  were  initially  chosen  based  on  those  for  bcc  tantalum  and  tantalum-tungsten  alloys14  and 
assigned  values  of  2/3,  1,  107  s"\  and  107  s-1  respectively.  The  normalized  activation  energy  (saturation),  goes , 
and  the  mechanical  threshold  stress  associated  with  “saturation”  of  the  evolving  structure,  i.e.,  dislocation 
population,  &eso  ,  are  calculated  using  a  Fisher-type  plot  approach.  Equation  (9)  is  re-expressed  in  the  form: 

lnGes  =  In  CTeso - In -  (20) 

ft b3  goes  e 


Figure  10  is  a  graphical  representation  of  equation  (20)  with  the  parameters  given  in  Table  HI. 

The  threshold  saturation  stress  at  a  given  temperature  and  strain-rate,  &es ,  can  be  estimated  by 
extrapolating  the  tangent  modulus  response  curve,  i.e.,  0  -  ag  response,  to  a  zero  strain-hardening  rate,  i.e.,  0  =  0; 
for  example,  see  Figure  1 1.  The  contributions  of  athermal  hardening  and  intrinsic  strengthening,  both  of  which  are 
independent  of  structure  evolution  (strain),  are  subtracted  from  the  overall  flow  stress  to  arrive  at  a  value  of  cte  : 


H  HO  lJSE 


(21) 


where  a  is  given  by  the  experimentally  determined  stress-strain  data.  Based  on  Figure  10,  deSO  and  go&s  are 
estimated  as  790  MPa  and  0.1 12  respectively. 
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Figure  10.  Modified  Arrhenius  Plot  Illustrating  the  Temperature  and  Strain-Rate 
Sensitivity  of  Structure  Evolution  Stress ,  ocs 

[Threshold  stress,  oEso ,  and  the  normalized  activation  energy  value,  goes,  are  determinable  from  this  plot.] 


Table  III.  Parameters  Used  to  Calculate  g0ts  and  6 


'eSO 


ln(  6e5  ) 

kT/nb3  ln(i£S0/z) 

OeS  [MPa] 

e  [s'1] 

T(K] 

[l  [MPa] 

6.3969 

0.02252 

600 

0.001 

77 

71,239 

6.2146 

0.05597 

500 

0.001 

188 

69,975 

6.4297 

0.03016 

620 

2,500 

298 

68,499 

5.9915 

0.07250 

400 

0.1 

298 

68,499 

5.8579 

0.09063 

350 

0.001 

298 

68,499 

6.2146 

0.05372 

500 

2,800 

473 

66,063 

5.9915 

0.07985 

400 

3,000 

673 

63,240 

5.7038 

0.10844 

300 

3,100 

873 

60,404 
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Figure  11.  Tangent  Modulus  Response,  9-o  Plot,  Determined  from  Large  Strain 

Torsion  Response 

[Superimposed  is  an  extrapolation  to  a  zero  strain-hardening  level  (flow  stress  saturation),  characterized  by  ocs.] 


The  saturation-rate  parameter,  a,  in  the  Voce-like  tank  strain-hardening  rule  was  determined  as  a  =  2  based 
on  previous  experience  with  other  alloys  systems  and  then  subsequently  optimized  through  an  iterative  process. 
Similarly,  the  stage  II  strain  hardening  rate,  0O,  is  given  as  by  an  empirical  relationship: 

60  =  5102.4  -  2.0758  T  (22) 

where  T  is  temperature  in  kelvins  and  0O  is  in  units  of  MPa/strain.  Unlike  the  basis  in  thermally  activated 
deformation  of  the  mechanical  threshold  stress,  a  and  the  tank  hardening  rule  are  completely  empirical  in  nature. 
Admittedly,  these  descriptions  of  the  structure  evolution  are  the  weakest  aspects  of  the  strength  model. 

Adiabatic  heating  effects,  relevant  under  dynamic  (high  strain-rate)  loading  conditions,  are  represented 
within  the  MTS  strength  model  framework  via  the  deformation  temperature  in  strengthening  contributions 
affected  by  thermal  activation.  For  HY-100  steel,  adiabatic  heating  is  presumed  to  occur  if  the  strain  rate  is  in 
excess  of  e  ~  500  s-1.  Temperature  increases  due  to  adiabatic  heating  are  expressed  by: 

AT-—  Ja(e)-Je  (23) 

pCp 
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where  ¥  can  be  analogized  with  the  fraction  of  work  converted  into  heat,  as  opposed  to  producing  strain 
hardening.  For  HY-100  steel,  *F  is  assumed  equal  to  0.95.  The  density  of  HY-100  steel  is  being  assumed 
equivalent  to  iron,  such  that  p  =  7.86  Mg/m3,  and  the  temperature-dependent  specific  heat  is  being  given  by  the 
empirical  relationship:11 

CP  =  0.09278  +  7.454M0-4  •  T  +  (24) 

rjy  JL 

in  units  of  MPa  m3/Mg/K. 
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DISCUSSION 


Constitutive  Behavior  of  HY-100  Steel 

The  MTS  framework,  in  general,  well  characterizes  the  deformation  response  of  HY-100  steel  near  the 
onset  of  yielding  (Figure  12).  Both  temperature  and  strain  rate  sensitivities  of  the  yield  strength  are  captured,  and 
their  inter-relationship  is  described  by  an  Arrhenius  expression,  Equation  (3).  A  single  Arrhenius  expression  is 
employed  to  characterize  thermally  activated  deformation  over  the  investigated  temperature/strain-rate  regime, 
cf.,  a,  and  Equation  (19).  No  dominant  thermally  activated  process,  however,  is  identified  as  affecting 
dislocation  motion.  Rather  an  “effective”  obstacle  is  prescribed,  encompassing  the  Peierls  effect  at  low 
(cryogenic)  temperatures  and  solute  strengthening  at  ambient  (298  K)  and  elevated  temperatures,  through  the 
parameters,  got.  Pi,  qi,  i0i ,  and  6, . 


Figure  12.  Comparison  between  MTS  Prediction  and  Uniaxial  Compression  Data 
for  HY-100  Steel  (Quasi-Static  Response)  (Page  1  of  2) 
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Figure  12.  Comparison  between  MTS  Prediction  and  Uniaxial  Compression  Data 
for  HY-100  Steel  (Quasi-Static  Response)  (Page  2  of  2) 
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A  marked  deviation  between  experiment  and  MTS  prediction  is,  however,  observed  for  HY-100  steel  for  a 
test  condition  at  77  K  and  10~3  s-1  (Figure  12).  The  relative  difference  between  yield  strength  determined  by 
experiment  and  MTS  prediction  is  =  24%.  MTS  calculation  underpredicts  yield  strength.  Although  the  nature  of 
this  disparity  is  not  currently  understood,  potential  explanations  may  include  (a)  the  presence  of  a  distinct 
deformation  mechanism  that  is  not  sufficiently  represented  by  the  current  single  “effective”  obstacle  description 
and  (b)  the  role  of  slip  mode,  i.e.,  pencil  glide  versus  restricted  glide,  on  yield  strength.15  Both  rationale  require 
modifying  the  current  three-term  MTS  strength  model  approach,  either  through  incorporating  a  temperature/rate- 
sensitive  strength  term,  as  for  tantalum14  and  zirconium  ,  or  amending  a  Taylor-like  factor.  The  latter  has  been 
considered  employing  polycrystalline  plasticity  calculations.15'17'18  In  short,  these  calculations  suggest  low- 
temperature  deformation,  if  dominated  by  restricted  glide,  i.e.,  { 1 10}<1 1 1>  slip,  can  result  in  ~  15%  larger  yield 
strengths  than  those  predicted  if  assuming  low-temperature  deformation  occurring  via  <lll>-pencil  glide  slip. 
Although  not  accounting  for  the  disparity  between  experiment  and  prediction  in  full,  considering  the  geometric 
nature  of  slip  improves  agreement  between  experiment  and  is  incorporated  into  the  thermally  activated 
deformation  framework. 

Strain-hardening  behavior  is  less  well  characterized  than  is  yield  strength  by  the  MTS  strength  model.  In 
particular,  the  compressive  stress-strain  response  illustrated  in  Figure  13  exhibits  increasing  deviation  between 
model  prediction  and  experiment  at  increasingly  larger  deformation  levels,  i.e.,  at  large  plastic  strains  under  room- 
temperature  deformation  conditions.  The  nature  of  this  discrepancy  arises  from  an  a  priori  assumption  (and 
enforcement)  of  flow  stress  saturation.  HY-100  steel,  however,  exhibits  strain-hardening  behavior  that  decreases 
from  the  initial  stage  II  value,  0O,  to  an  apparent  relatively  constant,  finite  strain-hardening  rate.  This  constant 
strain-hardening  rate  is  often  referred  to  as  linear  or  “stage  IV”  strain  hardening.  Similar  linear  hardening 
behavior  has  been  observed  in  other  steels19  and  aluminum  alloys20  and  is  of  potential  importance  in  loading 
scenarios  involving  large  strain  deformation. 


Figure  13.  Large  Strain  Compression  Data  and  MTS  Prediction  Comparison 
[Employing  (a)  modified  Voce  hardening  rule  imposing  flow  stress  saturation  and 
(b)  modified  Voce  hardening  rule  with  linear  strain-hardening  component] 
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Large  strain  uniaxial  compression  response  of  HY-100  steel  tends  to  indicate  a  linear  hardening  rate  of 
~  200  to  300  MPa/strain  is  achieved  at  flow  stresses  greater  than  960  MPa.  This  finite  degree  of  strain  hardening, 
or  lack  of  flow  stress  saturation,  at  large  deformation  levels  gives  rise  to  the  lack  of  agreement  between  model  and 
experiment,  and  also  indicates  the  current  Voce-like  tanh  hardening  rule  may  not  adequately  capture  structure 
evolution.  In  lieu  of  a  complex  dislocation-based  description  of  strain  hardening  encompassing,  for  example, 
cross-slip  and  dislocation  accumulation  and  annihilation  within  low  energy  dislocation  cell  structures,21'23  an 
empirical  description  of  hardening  analogous  to  that  of  the  Voce  rule  is  proposed. 

An  alternative  strain-hardening  formulation  based  on  the  Voce-like  expression  currently  utilized  within  the 
MTS  strength  model  framework  is  formulated  to  describe  the  linear  strain-hardening  behavior  of  HY-100  steel. 
Linear  hardening  is  captured  by  a  Voce-like  expression  by  assuming  the  initiation  of  stage  IV  hardening  marks 
the  termination  of  stage  III  in  the  absence  of  flow  stress  saturation.  In  other  words,  rather  than  stage  II 
transitioning  via  stage  HI  to  flow  stress  saturation,  as  currently  described  in  the  MTS  strength  model,  stage  H 
transitions  via  stage  HI  to  linear,  stage  IV  hardening.  A  modification  to  an  existing  tanh  strain  hardening  rule 
(equation  (8)  is  proposed  to  provide  (a)  a  general  method  to  represent  either  flow  stress  saturation  or  stage  IV 
(linear  hardening)  behavior,  (b)  an  improved  stress-strain  prediction  at  both  small  and  large  strain,  and  (c)  a 
formulation  that  can  be  readily  integrated  to  yield  a  closed-form,  analytic  solution  of  e  =  fia).  The  linear 
hardening  rule  is  imposed  such  that  achieved  is  a  continuous  increase  in  flow  stress  with  increasing  strain  and  a 
constant  hardening  rate  at  large  strains  (and  flow  stresses)  (Figure  13). 

A  proposed  “1st  order”  modification  to  the  tanh  strain-hardening  rate  expression  (with  linear  stage  IV 
strain-hardening)  is  given  below  as: 
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Analogous  to  equation  (8),  the  modified  strain-hardening  rate  term  (equation  (25)  is  composed  of  hyperbolic 
tangent  terms.  This  modified  hardening  expression  takes  advantage  of  the  former  tanh  rule’s  good  description  of 
stress-strain  response  at  small  strains  levels.  Furthermore  stage  IV  hardening  at  large  strain  levels  can  be  better 
represented  by  the  linear  hardening  rate  modification. 

The  contribution  due  to  linear  (stage  IV)  hardening  is  present  as  the  second  term  on  the  right  hand  side  of 
equation  (25).  Saturation  stress,  as  present  in  the  Voce-like  (e.g.,  tank)  expressions,  is  effectively  redefined  as  the 
stress  level  at  which  flow  stress  saturation  would  occur  in  the  absence  of  linear  hardening.  In  the  absence  of  stage 
IV  hardening,  the  saturation  stress,  as  originally  intended,  in  the  Voce-like  (e.g.,  tank)  expression  can  be 
introduced  into  equation  (25)  and  the  value  for  the  stage  IV  strain-hardening  rate,  0iV,  set  equal  to  zero. 
Conceptually,  the  above  expression  can  be  envisaged  as  stage  HI  behavior,  as  represented  by  the  first  tanh(...) 
term,  provides  a  transition  between  a  first  regime  of  linear  strain-hardening,  i.e.,  stage  II  (0O),  to  a  second  regime 
of  much  smaller  strain-hardening,  i.e.,  stage  IV  (0|V). 


Identical  to  equation  (8)  equation  (25)  can  be  integrated  directly  with  respect  to  6e  to  yield  to  a  closed- 
form  analytic  solution  for  strain: 
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which  can  assessed  over  a  finite  stress  interval  between  ctg  j  and  ctg  2  •  Recall  again  that  within  the  MTS 
framework  stress  is  the  internal  state  variable  and,  as  such,  e  =j^0).  As  mentioned  previously  dcs  is  redefined  to 
characterize  the  threshold  stress  at  which  stage  III  hardening  transitions  to  stage  IV,  in  the  absence  of  flow  stress 
saturation.  Parameters  for  this  modified  hardening  rule  (equation  26)  are  shown  below: 
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Torsion  Versus  Uniaxial  Compression  and  Texture  Evolution 

Uniaxial  compression  and  torsion  tests  provide  constitutive  data  from  relatively  well-defined  stress-states, 
e.g.,  uniaxial  or  pure  shear,  respectively,  and  are  both  utilized  when  available.  Figure  14  illustrates  that  the 
mechanical  responses  derived  from  uniaxial  compression  and  torsion  are  in  good  agreement  up  to  an  effective 
plastic  strain  of  -0.40.  However,  increasing  deviation  in  effective  stress  values  are  observed  with  increasing 
plastic  strain.  The  compressive  stress-strain  response  appears  to  exhibit  near-linear  hardening  at  large  strain 
levels,  whereas  the  torsion-derived  stress-strain  behavior  appears  to  approach  flow  stress  saturation. 


Figure  14.  Comparison  between  Large  Strain  Compression  and  Torsion  Data 
[Increasing  deviation  between  response  is  observed  with  increasing  deformation.] 


Polycrystalline  plasticity  calculations,  using  the  LApp  code,6  suggest  deformation-induced  texture  evolution 
results  in  a  “harder”  material  when  the  deformation  path  is  uniaxial  compression,  as  opposed  to  torsion.  ’ 
Figures  15  and  16  illustrate  the  Taylor  and  texture  factor  evolutions  as  a  function  of  deformation  for  a  uniaxial 
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compression  and  torsion.  The  increase  in  the  Taylor  factor  from  -2.73  to  -2.85  (e  -  1.0)  in  compression  indicates 
textural  hardening  occurs.  In  pure  shear  a  decreasing  Taylor  factor  from  -2.73  to  -2.66  (s  -  1.0)  is  observed  and 
tends  to  indicate  textural  softening  occurs.  Furthermore,  the  calculation  results  illustrate  that  differences  in  Taylor 
factor  between  compression  and  shear  are  exacerbated  with  increasing  deformation;  e.g.,  MCOmpression/Mt0rsion 
increases  with  increasing  deformation.  These  contrasting  Taylor  factor  evolutions  in  part  rationalize  the 
differences  in  stress-strain  response  observed  between  uniaxial  compression  and  torsion  with  increasing 
deformation. 

Stress-strain  curves  (Figure  17)  illustrate  that  compensating  for  deformation-induced  texture  evolution 
brings  into  better  agreement  the  uniaxial  compression  and  torsion  mechanical  responses;  the  procedure  used  to 
account  for  texture  evolution  is  presented  in  Appendix  B.  Qualitatively,  the  similar  strain-hardening  rates 
exhibited  in  compression  and  torsion  suggest  the  hardening  mechanisms  are  similar  between  loading  paths  and 
can  be  modeled  with  an  identical  hardening  rule  and  parameters  by  the  MTS  strength  model. 

The  linear  strain-hardening  modified  MTS  strength  model  prediction  of  the  room  temperature,  quasi-static 
compression  data  is  re-evaluated  to  include  texture  evolution  effects.  Taylor  factor  evolution  is  incorporated  into 
the  structure  evolution  term,  such  that: 


= 


* 


(27) 


where  M*  is  the  normalized  Taylor  factor,  as  given  in  Appendix  B.  MTS  parameters  correspond  to  those 
presented  above.  These  MTS  strength  model  parameters,  as  well  as  those  listed  in  Appendix  A,  can  be  considered 
as  texture  compensated  and  may  be  applied  in  representing  effective  stress-strain  response  in  conjunction  with 
appropriate  texture,  i.e.,  Taylor  factor,  evolution  rules,  e.g.,  equations  (B2)  and  (B3),  for  the  current  HY-100  steel 
plate  material.  Figure  18  illustrates  MTS  strength  model  predictions  with  and  without  texture  evolution  and  the 
effect  of  texture  evolution  (softening)  under  a  torsional  loading  path.  This  figure  also  shows  that  incorporating 
texture  evolution  into  the  constitutive  model  may  explain  in  part  the  apparently  dissimilar  stress-strain  curves 
measured  in  uniaxial  compression  (from  which  the  MTS  parameters  are  determined)  and  torsion.  However,  the 
fair  agreement  between  torsion  and  the  MTS  prediction  illustrates  accurate  descriptions  of  texture  evolution  in 
torsion  and  compression  are  still  required  as  well  as  a  more  robust  implementation  of  texture  evolution  in  the 
strength  model. 
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Figure  15.  Taylor  Factor  Evolution  under  Uniaxial 
Compression  and  Torsional  Loading  Paths 
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Figure  16.  Texture  Evolution  Manifested  by  Inverse  Pole  Figures  as  a  Function  of 

Deformation  and  Loading  Path 
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Figure  17.  Stress-Strain  Curves  Normalized  for  Loading-Path-Dependent 

Taylor  Factor  Evolution 


Figure  18.  Compensating  for  Texture  Evolution  in  the  MTS  Strength  Model 
to  Rationalize  the  Dissimilar  Stress-Strain  Responses  Obtained  in  Compression 

and  Torsion 
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CONCLUSIONS 


The  MTS  strength  model  framework  is  well-suited  to  characterizing  the  constitutive  mechanical  response  of 
metals  and  alloys.  A  three-term  MTS  strength  model  currently  characterizes  Navy  HY-100  steel.  Featured  within 
the  three-term  MTS  model  are  constant  structure  considerations,  such  as  athermal  and  thermally  activated 
deformation  mechanisms,  and  structure  evolution  (strain-hardening)  considerations. 

Athermal  strengthening  is  presumed  associated  with  martensite/bainite  lath  boundary  or  Hall-Petch-type 
behavior.  Thermally  activated  processes  encompass  the  intrinsic  Peierls  barrier  and  solute  strengthening  due 
principally  to  interstitial  atoms,  such  as  carbon,  nitrogen,  and  oxygen.  Additionally,  “constant  structure”  and 
thermally  activated  deformation  includes  the  initial  dislocation  population  within  the  tempered  martensite 
laths.  Structure  evolution  involves  the  dynamic  processes  of  dislocation  accumulation  and  annihilation. 
Analogous  to  the  concept  of  thermally  activated  deformation  of  the  “constant  structure,”  temperature/rate  effects 
are  modeled  for  structure  evolution  via  an  Arrhenius  expression.  However,  structure  evolution,  viz.,  the 
competition  between  dislocation  accumulation  and  annihilation,  is  modeled  solely  from  an  empirical  basis  and 
utilizes  a  Voce-like  expression.  The  Voce-like  strain-hardening  (structure-evolution)  expression  appears  to  be 
deficient  in  accurately  modeling  material  mechanical  response  at  large  deformation  levels  and  as  such  is  a 
principal  limitation  of  the  MTS  strength  model. 

A  principal  motivation  to  modify  the  existing  tank  (Voce-like)  hardening  rule,  as  currently  implemented  in 
the  MTS  strength  model,  is  to  improve  the  predicted  constitutive  stress-strain  behavior  at  both  small  and  large 
strain  levels,  with  primary  emphasis  on  improving  the  large  strain  stress-strain  prediction.  The  existing  tank 
hardening  rule  well  approximates  a  material’s  small  strain  ( s  <  0.10)  stress-strain  response.  However,  the  fidelity 
of  this  tank  rule  is  degraded  at  larger  strain  levels  (e  >  0.20  -  0.30).  A  modification  to  this  hardening  rule  is 
proposed  to  improve  the  large  strain  prediction  while  retaining  the  good  small-strain  prediction  characteristics  of 
the  original  tank  model. 

Agreement  among  effective  stress-strain  data  obtained  from  different  loading  paths  is  improved  when 
experimental  stress-strain  data  is  compensated  for  deformation-induced  texture.  The  current  MTS  framework, 
equations  (1)  through  (11),  does  not  include  texture  considerations  in  the  structure  evolution  rule,  i.e.,  the  tank 
model.  As  such,  the  MTS  strength  model  incorporates  both  dislocation  hardening/softening  and  texture  evolution 
into  the  structure  evolution  rule.  Polycrystalline  plasticity  calculations  illustrate  that  deformation-induced  texture 
evolution,  manifested  by  a  strain-dependent  Taylor  factor,  becomes  marked  at  deformation  levels  of  8  >  0.20. 
Furthermore,  notable  deviation  between  effective  stress-strain  responses  is  observed  between  uniaxial 
compression  and  torsion  in  the  absence  of  compensating  deformation  for  texture  evolution.  Subsequent  to 
considering  texture  evolution  in  the  stress-strain  response,  very  good  agreement  is  achieved  between  uniaxial 
compression  and  torsion. 

An  implication  of  the  current  observations  and  results  is  that  improved  agreement  between  effective  stress- 
strain  response  over  arbitrary  loading  paths  can  be  achieved  if  texture  evolution  is  compensated  for  within  the 
experimental  data  and  strength  model.  The  example  presented  above  illustrated  the  agreement  between  effective 
stress-strain  behavior  for  a  near-random  HY-100  steel  deforming  over  a  well-defined  loading  path  involving 
either  uniaxial  compression  or  torsion  when  texture  evolution  was  coupled  to  structure  evolution.  Highly  textured 
materials,  such  as  tantalum,  may  exhibit  significant  deformation-induced  texture  evolution  and  may  require  that 
for  accurate  stress-strain  representation  the  constitutive  strength  model  include  a  structure  evolution  term  to 
account  for  Taylor  factor  evolution. 
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Appendix  A 


MTS  EQUATIONS  AND  PARAMETERS  FOR  HY-100  STEEL 
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Where:  po  =  7. 146  x  104  MPa 
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b 3 

it  =  1.38  x  10-23  j/k,  b  =  2.48  x  10~10 
E0i  =  1  x  1013  s_1 

8oi=  1161 
q{  =  3/2,  Pi  =  1/2 

Where:  e0E  =  lx  107  s_1 

8oe  =  16 
*?£  =  1>  Pe  “  2/3 


Where: 

00  =5102.4-  2.0758  -T;MPa 
a  =  2 


\^_Soes_ jn  ®es 

kT 


Where:  0.112 

o  csO  —  7  82  MPa 
eoej  =  lx  107  s"1 
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.  /IT 


*p 

-—fa(e)de; 

pCp 


when  £>  500  j  1 


Where:  ¥  =  0.95 
p  =  7.86  Mg/m3 

Cp  =  0. 09278  +  7. 454  x  10~4  •  T  + 
MPa  m3/Mg/K 


12404 

T2 


A-4 


Appendix  B 

TEXTURE  EVOLUTION  COMPENSATION 
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Texture  evolution  in  the  experimental  data  is  accommodated  by  normalizing  the  experimentally  determined 
flow  stress  by  a  deformation-dependent  relative  Taylor  factor.  The  differences  in  stress-strain  response  amongst 
loading  paths  are  addressed  by  compensating  experimental  uniaxial  compression  and  torsion  data  for  texture,  cf., 
equation  (17),  and  texture  evolution.  Experimental  data  are  normalized  for  texture  evolution  through: 


^  NoTcxture  _  q 


(Bl) 


where  CTNoTexture  corresponds  to  texture-compensated  experimental  data,  aExpt  refers  to  the  original  experimental 
data,  and  M*  is  the  appropriate  normalized  Taylor  factor  and  equal  to  M(e)/M0.  The  M-evolution  over  a  uniaxial 
compression  loading  path  is  expressed  in  terms  of  a  sixth-order  polynomial  (cf.,  Figures  15  and  16): 

M(s)compress i0„  =  -1.3 14 le6  +  3.6953e5  -  3.2263s4  +  0.369s3  +  0.6234e2  -  0.0297s  +  2.7302  (B2) 

Similarly,  the  evolution  of  M  over  a  pure  shear  loading  path  is  expressed  as  (cf.,  Figures  15  and  16): 

M(e)torsion  =  -0.1632s6  -  0.5208s5  +  2.1803s4  -  2.0809s3  +  0.5296s2  -  0.0184s  +  2.7306  (B3) 

The  initial  Taylor  factor  at  s  =  0  is  M0  =  2.73.  The  resulting  CTNoTexture  -  e  curves  are  illustrated  in  Figure  16. 
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